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We present a calculation of the renormalisation constants for the temporal vector current, Zy and spatial axial 
current, ZA k , to O(jf) for B — > D transitions using the O(^) NRQCD action for both fo and c quarks evaluated 
for a large range of mass parameters. Considerations for the renormalisation of the spatial vector current and the 
temporal axial current are discussed and initial results for a mixed lattice current are presented for the spatial 
vector current. 



1. Introduction 

The semi-leptonic B — ► D, D* decay is phe- 
nomenologically interesting since its Feynman 
amplitude involves the poorly know CKM matrix 
element V c b 

M{B - XI- v) = - % -V ch u a »{l - ls)v v Hp (1) 

Hp = (X(p')\c^(l- l5 )b\B(p)) , X = D,D*(2) 

where the hadronic tensor H^, and its Lorentz 
decompostion into form factors is calculable using 
lattice QCD. 

We perform the matching of the Lattice to con- 
tinuum MS currents to 1 loop using the on-shcll 
scheme in Feynman gauge with vanishing exter- 
nal spatial momenta. In the on-shell scheme the 
wavefunction renormalisation and vertex correc- 
tions individually contain (cancelling) infra-red 
divergences which we control using a fictitious 
gluon mass, A. We use dimensional regularisa- 
tion to control the ultra-violet behaviour in the 
continuum calculation. 
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Here the currents Vq and A k can be renor- 
malised to O(m) m ^ ne usua l manner, however 
the currents V& and Aq will involve operator mix- 
ing at this order. 

3. 1-loop Lattice Correction 

We use the 0(±) lattice NRQCD action § for 
both b and c quarks, 
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2. Continuum Calculation 



Expanding the 1-loop correction to the currents 
to O(-lj), and writing e 



xr-, we obtain 
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for which the Feynman rules may be found in 
reference 

The Foldy-Wouthysen transformation of the 
continuum current corrections and dimensional 
arguments for which currents can contribute at 
0(jj) suggest we take the bases of operators in 
Table |l| for the lattice currents. 

For those currents that are unmixed, the calcu- 
lation is similar to , and we write the renormal- 
isation constant Zj? = 1 



a s Zp = 1 
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Table 1 

Lattice Current Bases 
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computed numerically using VEGAS ||, per- 
forming the temporal loop momentum integra- 
tion analytically for those contributions that were 
infra-red divergent. 

In the case of the spatial vector and tempo- 
ral axial currents it is necessary to evaluate the 
derivatives of the diagrams with respect to the ex- 
ternal momentum to match with the continuum. 
This is done by numerically evaluating the inte- 
gral of the analytically taken derivative, resolving 
the pieces of the Pauli structure by taking differ- 
ent derivative directions. 

In the continuum, both the vertex correction to 
the current cTb, and the wavefunction renormali- 
sation contained a logarithmic divergence propor- 
tional to log -pT where the sign is positive in 
the vertex correction. The same infra-red diver- 
gences were found in the lattice vertex correction 
and wavefunction renormalisation so that as one 
would expect the infra-red behaviour of the the- 
ory is identical to that in the continuum. 

3.1. Vq and A k Currents 

The 1-loop contribution to the renormalisation 
constant Zy is plotted in Figure ^ for various 
values of e = -jgjj- and M&. Here the stabilisation 
parameter n has been chosen in a mass dependent 
manner such that Mo xn > 3. The lattice current 
is the conserved current of the lattice action, and 
the continuum current is conserved, so that there 
is no correction on the line e = 1. 

The 1-loop correction to the spatial axial cur- 
rent is in Figure ||. The mass dependence is 
stronger since the line for e = 1 is not protected. 

3.2. 14 Current 

We define the 1-loop lattice mixing matrix Zij , 



(c(p')\j i \b( P )} = j2^. 
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Figure 1. Z$ as a function of -tj- and ¥^ L . The 



lines correspond to If = 0.1,..., 1.0 in incre- 



ments of 0.1, with the = 0.1 curve topmost. 

where the JT^ are the continuum analogs of the 
lattice operators given in Table g. We define the 
contribution to the mixing matrix arising from 
the lattice vertex correction, fy, via 
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We then invert the mixing matrix and match to 
the continuum. Here the coefficients Bi may be 
inferred from equation [|, and the Z m factors arise 
from matching the bare mass in lattice currents to 
the pole mass in the tree level continuum current. 
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The diagonal elements £jj contain infra-red di- 
vergences which cancel with those in the lattice 
wavefunction renormalisations. Those currents 
appearing at tree level also carry tadpole im- 
provement counter terms. We therefore rewrite 
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Spatial Axial Current 
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Figure 2. Z 1 } 1 as a function of ^ and The 
lines correspond to -jgjj- = 0.1,. ...1.0 in incre- 
ments of 0.1, with the j^- = 0.1 curve topmost. 

Equation [ll] in terms of IR divergence and tad- 
pole counter term free quantities, as follows: 
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Preliminary calculations of the mixed lattice 
have been performed on a few mass values for 
degenerate quarks. The degenerate case may in 
fact be of use in certain lattice simulations to 
obtain the Isgur-Wise function, however the re- 
sults should really be considered illustrative of the 
method, and the calculation over a similar param- 
eter regime to our previous results for the other 
currents will be performed. In this case the coeffi- 
cients of Ji and J2, and of J3 and J4 are identical, 
and we only present the coefficients for J\ and J3 

in Table 0. We evaluate the tadpole improvement 

[11 

counter terms using u (plaq) = -| , and take val- 
ues for Z m from Also, for the degenerate case 
Bi = B 2 = —J- and = Ba = —J^, so that 

1 Z 37T ° ^ 37T ' 



Table 2 

One Loop Correction Coefficients 
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£13 + £23 


Zip 


2 


-0.6581(5) 


-0.2516(1) 


-0.4156(6) 
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-0.9219(5) 


-0.2041(1) 


-0.0030(6) 


10 


-1.2662(5) 


-0.0963(1) 


-0.2711(8) 
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One Loop Correction Coefficients 
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2,2 


-0.776(4) 


0.145(1) 




4,2 


-0.571(4) 


0.098(1) 




10,1 


-0.247(4) 


-0.010(1) 





when the mixed current is written as 

Vfs = (l + a,Ci)(J 1 + Ja) 

+ a s C 2 (J 3 + J 4 ), (13) 

we obtain the coefficients in Table [| 

4. Future Direction 

We are nearing completion of the calculation 
of the mixed current for 14, while there is a 
smaller basis of operators for Aq, so that the cor- 
responding calculation should be somewhat eas- 
ier. Thereafter we shall perform a similar calcu- 
lation for the b — > c transition using the NRQCD 
action for the b quark, and the O(a) improved 
Wilson action for the c quark. The extension of 
the calculation to 0(jj?) would be of interest. 
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